home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DGBSL.z / DGBSL
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDGGGGBBBBSSSSLLLL((((3333FFFF))))                                                            DDDDGGGGBBBBSSSSLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DGBSL   - DGBSL solves the double precision band system A * X = B  or
  10.      TRANS(A) * X = B using the factors computed by DGBCO or DGBFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE DGBSL(ABD,LDA,N,ML,MU,IPVT,B,JOB)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAABBBBDDDD DOUBLE PRECISION(LDA, N)
  20.         the output from DGBCO or DGBFA.
  21.  
  22.      LLLLDDDDAAAA INTEGER
  23.         the leading dimension of the array  ABD .
  24.  
  25.      NNNN INTEGER
  26.         the order of the original matrix.
  27.  
  28.      MMMMLLLL INTEGER
  29.         number of diagonals below the main diagonal.
  30.  
  31.      MMMMUUUU INTEGER
  32.         number of diagonals above the main diagonal.
  33.  
  34.      IIIIPPPPVVVVTTTT INTEGER(N)
  35.         the pivot vector from DGBCO or DGBFA.
  36.  
  37.      BBBB DOUBLE PRECISION(N)
  38.         the right hand side vector.
  39.  
  40.      JJJJOOOOBBBB INTEGER
  41.         = 0         to solve  A*X = B ,
  42.         = nonzero   to solve  TRANS(A)*X = B , where
  43.         TRANS(A)  is the transpose.  On Return
  44.  
  45.      BBBB the solution vector  X . Error Condition
  46.  
  47.      AAAA division by zero will occur if the input factor contains a zero on the
  48.      diagonal.  Technically this indicates singularity but it is often caused
  49.      by improper arguments or improper setting of LDA .  It will not occur if
  50.      the subroutines are called correctly and if DGBCO has set RCOND .GT. 0.0
  51.      or DGBFA has set INFO .EQ. 0 .  To compute  INVERSE(A) * C  where  C  is
  52.      a matrix with  P  columns
  53.         CALL DGBCO(ABD,LDA,N,ML,MU,IPVT,RCOND,Z)
  54.         IF (RCOND is too small) GO TO ...
  55.         DO 10 J = 1, P
  56.         CALL DGBSL(ABD,LDA,N,ML,MU,IPVT,C(1,J),0) 10 CONTINUE LINPACK.  This
  57.      version dated 08/14/78 .  Cleve Moler, University of New Mexico, Argonne
  58.      National Lab.  Subroutines and Functions BLAS DAXPY,DDOT Fortran MIN0
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.